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ABSTRACT 



The Milky Way has a giant stellar structure in the solar neighborhood, which 
has a size of ~ 1 kpc, a mass of ~ IO^Mq, and a ring-like distribution of young stars. 
Fundamental physical properties of this local enigmatic structure, known as the Gould 
belt (GB), have not been reproduced by previously proposed models. We first show 
that the local enigmatic structure can be formed about 30 Myr ago as a result of a 
high-speed, oblique collision between a gas cloud with a mass of ^ IO^Mq and a dark 
matter clump with a mass of ~ 1O''M0 based on numerical simulations of the collision. 
We find that strong dynamical impact of the clump transforms the flattened cloud 
into a ring-like stellar structure after induced star formation within the cloud. Our 
simulations furthermore demonstrate that the stellar structure is moderately elongated 
and significantly inclined with respect to the disk of the Milky Way owing to the 
strong tidal torque by the colliding clump. We thus suggest that the GB is one of 
stellar substructures formed from collisions between gas clouds and dark matter clumps 
predicted in the hierarchical clustering scenario of galaxy formation. We also suggest 
that collisions of dark matter clumps with their host galaxies can significantly change 
star formation histories for some of their gas clouds thus influence galactic global star 
formation histories to some extent. Our simulations predict that unique giant stellar 
substructures recently discovered in other galaxies can result from dynamical impact 
of their dark matter clumps on their gas clouds. 

Key words: The Galaxy - galaxies:structure - galaxies:kincmatics and dynamics - 
galaxies :halos - galaxies:star formation 



1 INTRODUCTION in the solar neighborhood (e.g., Poppel 1997). Other studies 

r, , (e.g., Comeon & Torra 1994) proposed an oblique impact of 

The GB is a flattened stellar structure with an inclination , . , , , , . ,, ^ , ,. , r ,i 

r 1 „„o , ^ , . , n . , ^^ny ^ high- velocity cloud on the gas m the Galactic disk tor the 

oi 16 — 22 to the Galactic plane and contains about oOyo r~.n r 

GB tormation. 

OI young stars with ages oi 30 — 60 Myr withm 600 pc irom 

the Sun (e.g., Olano 1982; Poppel 1997; GuiUout et al. 1998; The purpose of this paper is thus to propose a new sce- 

Torra et al. 2000). The projected distribution of young stars nario which explains self-consistently the observed funda- 

in the GB is suggested to be quite inhomogeneous and sig- mental properties of the GB. In the new scenario presented 

nificantly elongated with an eUipticity of ~ 0.3 (e.g., Poppel in this paper, the GB is one of stellar substructures (e.g., 

1997). The observed stellar kinematic of the GB suggests "ngs and arcs) formed from high-speed, off-center collisions 

that the GB has internal rotation with a possible angular between giant molecular clouds, which are believed to be the 

velocity of 7 - 24 km s"^ kpc'^ and is currently expanding birth places for stars (Lada & Lada 2003), and dark matter 

(Lindblad 2000). The GB is considered to be young with its clumps (DMCs) orbiting the Galaxy. This dynamical impact 

age less than 30 - 60 Myr (e.g., Poppel 1997). Previous the- of DMCs is referred to as the dark impact for convenience 

oretical studies suggested that the formation of the GB is in Bekki & Chiba (2006) and used as such in the present 

closely associated with eflicient star formation in an expand- Paper. Since recent high-resolution cosmological simulations 

ing gaseous ring formed as a resuh of stellar feedback effects of the Galactic halo formation based on the cold dark mat- 

of massive OB stars and supernovae on interstellar medium ter model (CDM) have shown that hundreds of high-density 

DMCs may exist in the solar-neighborhood (Diemand et al. 

2008), we consider that such collisions between gas clouds 

* E-mail: bekki@phys.unsw.edu.au and DMCs are almost inevitable. We here do not intend to 





Figure 1. Morphological evolution of a gas cloud colliding with a dark matter clump (DMC) projected onto the X-Y plane (upper) 
and the X-Z one (lower). The time T (in Myr) shown in the upper left corner of each panel represents the time that has elapsed since 
the simulation started. The two separate panels show the orbits of the cloud (cyan) and the DMC (green) projected onto the X-Y and 
the X-Z planes. The final locations of the cloud and the DMC indicated by an open square and and an open triangle, respectively, and 
the solar radius (_R = 8.5 kpc corresponding to the distance between the Sun and the Galactic center, G.C.) and the Galactic plane are 
shown in these two panels for comparison. The cloud and the DMC are represented initially by 150000 gaseous particles and 100000 
stellar ones, respectively, and shown by magenta and blue, respectively. The yellow particles are new stars formed from gas in the cloud. 
A bar shown in the lower left corner of each frame measures 100 pc. The time evolution of the cloud-DMC distance is described in Fig. 
3, which clearly shows that the cloud-DMC collision is an off-center one. 



discuss the well known over-abundance problem of satellite 
galajcies (hosted by dark matter subhalos) in the CDM, just 
because this is beyond the scope of this paper. 



2 THE MODEL 

We here show the results of GRAPE-SPH simulations (Bekki 
& Chiba 2007) only for the best model explaining the ob- 
served properties of the GB self-consistently, and will de- 
scribe those of other models with different model parameters 



in our future study (e.g., Bekki 2009 in preparation, B09). 
The progenitor cloud of the GB is modeled as an oblate 
sphere with a radius of 200 pc, an ellipticity of 0.3, a mass of 
IO^Mq, temperature of lOOK, and rigid rotation with an an- 
gular velocity of 4.8 x 10~^ km s~^ pc~^. These parameters 
values are consistent with the observed properties of GMCs 
(e.g., Solomon et al. 1979; Kerton et al. 2003; Rosolowsky et 
al. 2003). GMCs with sizes of ~ lOOpc exist in the Galaxy 
(e.g.. Fig. 6 in Solomon et al. 1979), and those with sizes 
of 100 — 200pc exist in other galaxies (e.g., NGC 4449 and 
NGC4605; see Fig. 1 in Bolatto et al. 2008). Therefore, given 
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Figure 2. Final distribution of stellar particles in the cloud at T = 45 Myr in the Galactic coordinate system, where abscissa and ordinate 
denote Galactic longitude (Z) and latitude (6), respectively. All stellar and gaseous particles are shown by red and the final location of the 
DMC is shown by a filled blue circle for comparison. The five blue open circles represent the observed OB associations on the Gould Belt 
(de Zeeuw et al. 1999): two OB2 in Scorpius constellation (Sco), OB 1 in Orion (Ori), 0B2 in Perseus (Per), and 0B2 in Cepheus (Cep). 
Only one of the listed locations in de Zeeuw et al. (1999) for each OB association is shown here for convenience (e.g., (i,fe)=(164°, — 13°) 
for Per OB2). The dotted thick cyan line describes the best-fit Gould curve (Taylor et al. 1987) of fe = 6osin( — ^o]) + ''c for the 
simulated stellar belt, whore bo (^ 12°), Iq (~ 232°), and 6c {~ —2°) denote the inclination angle of the belt with respect to the Galactic 
plane, the ascending node longitude of the belt, and the vertical distance between the mass center of the belt and the Galactic plane, 
respectively. The simulated Gould curve is more similar to that derived from the observed molecular clouds (Taylor et al. 1987) than to 
those from young stars (Poppel 1997; Perrot & Grenier 2003). 



the observed uncertainties in estimating the GMC sizes (Bo- 
latto et al. 2008), the adopted cloud sizes are consistent with 
observations for very large GMCs. Star formation is assumed 
to happen in local regions of the cloud if the regions satisfy 
physical conditions for the Jeans instability that can cause 
collapse of interstellar gas and subsequent star formation 
(e.g., Bekki & Chiba 2007). 

The self-gravitating cloud is dynamically influenced by 
the fixed gravitational potential of the Galaxy composed of 
bulge, disk, and dark matter halo. We adopt the same Galac- 
tic potential as that used in Helmi & de Zeeuw (2000) which 
investigated orbits of dwarfs galaxies in the Galactic halo. 
The cloud is located initially at 8.5 kpc from the Galactic 
center (G.C.) and moves round the Galaxy with a circular 
speed of 220 km s~^. The present location of the Sun is set 
to be {X, Y, Z) — {~8.5, 0, 0) kpc and the Galaxy rotates 
to the +Y direction in the present study. The initial loca- 
tion and velocity for the cloud are chosen such that the final 
(i.e., present) location can be almost the same as that of 
the Sun after its 45 Myr evolution (see Fig.l for the initial 
three-dimensional position in the Galaxy). 

The DMC has a mass of 3 X IO^'Mq, a mean density 
of 2.OM0 pc"^ within lOOpc, and the NFW radial density 
profile (Navarro et al. 1996) within its tidal radius {— 317 
pc), which are all consistent with predictions from the lat- 
est cosmological simulations (Diemand et al. 2008; Springel 
et al. 2008a). The initial relative position {Xi, Yi, Zi) and 
velocity {Ur, Vi, Wt) of the DMC with respect to the cloud 
are set to be (-446, 233, -600) pc and (59, -29, 66) km 
s~^, respectively. For these relative positions and velocities. 



the DMC collides with the cloud ai T ^ 8 Myr, where T 
represents the time that has elapsed since the simulation 
started. 

We have confirmed that (i) the ring-like structure can 
be formed in models with a wider range of model param- 
eters (e.g., impact parameter) and (ii) the dark impact of 
DMCs with masses of ~ IO'^Mq can not influence struc- 
ture and kinematics of the old stellar disk (B09). The origin 
of unique giant stellar substructures like the GB found in 
disks for a number of galaxies (Efremov& Elmegreen 1998; 
Comeron 2001; Larsen et al. 2002) and the "shingles" in the 
Galaxy (e.g., Schmidt-Kaler & Schlosser 1973 ) can be also 
explained in the context of the dark impact (B09) . The sim- 
ulated inclination angle of the GB is slightly smaller than 
the observed one, which may be one of disadvantages of the 
model: however, the inclination angle depends on (i) the 
time step at which it is estimated and (ii) the assumed lo- 
cation of the Sun owing to the presession of the GB. 



3 RESULTS 

Fig. 1 shows that as the DMC passes through the cloud 
from its beneath (T = 7 Myr), it first compresses the cen- 
tral region of the cloud so that new stars can form from 
the high-density region of the cloud (T = 7 — 15 Myr). Star 
formation rate increases rapidly during the collision (T = 8 
Myr), reaches its maximum of 0.16 M© yr"^ (T = 12 Myr), 
and rapidly decrease to be 0.005 Mq yr~^ after the collision 
(r = 22 Myr). The dark impact then induces an outwardly 
propagating density wave and consequently transforms the 
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Figure 3. Time evolution of the distance between the gas cloud 
and the DMC (left) and the projected distribution of new stars 
in the simulated stellar structure at T = 45 Myr (right). The 
dotted thick green line indicates the initial cloud radius (=200 
pc) and the intersection between the line and the thick solid green 
line for the cloud-DMC distance can indicate the time when the 
cloud and the DMC start their strong dynamical interaction in 
the left panel. The dotted red lines indicate the start and the end 
of active star formation (i.e., starburst) during the dark impact 
in the left panel. The two dotted blue lines indicate a range of 
younger stellar populations with ages (tage) between 1 Myr and 
10 Myr in the left panel. The starburst population with 25 Myr 
^ iage ^ 37 Myr and younger one with 1 Myr ^ tage ^ 10 Myr are 
shown by smaller red dots and bigger blue ones, respectively in 
the right panel. Here AX and Ay denote X- and K-positions of 
particles with respect to the mass center of the stellar structure. 
The direction of the Galactic center (G.C.) is shown by a green 
arrow. The center of the frame is set to be the mass center of 
the stellar structure so that the simulated distribution can be 
compared with the observed one (Poppel 1997; Perrot &; Grenier 
2003). 
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Figure 4. Kinematics of the simulated stellar structure. The 
radial velocities (vi) and X-components of velocities (■Ux) for new 
stars with 25 Myr ^ tago ^ 37 Myr (smaller red dots) and with 
1 Myr ^ tagc ^5 10 Myr (bigger blue dots) are plotted against 
Ay in upper and lower frames, respectively. Here and Dx are 
velocity components with respect to those of the mass center of 
the stellar structure. The dotted green lines in the upper and 
lower panels represent the mean value of Vr (= 7.1 km s~^) for 
all stellar particles and rigid rotation with an angular velocity of 
21.1 km s~^ kpc~^, respectively. 



oblate gas cloud into the ring-like one (T — 22 Myr). The 
tidal torque of the dark impact also causes the precession 
of the internal spin vector of the new stellar system so that 
the system can be seen as significantly inclined with respect 
to the Galactic plane (T — 22 Myr). The gas and new stars 
in the cloud can continue to expand after the dark impact, 
while the inclination angle of the stellar system becomes ap- 
parently smaller in edge-on view owing to precession {T — 30 
Myr) . The strong tidal field of the Galaxy transforms the ex- 
panding stellar system into an elongated ring-like structure 
that is still inclined to the Galactic plane (T=45 Myr). 

Fig. 2 shows that the simulated ring-like structure can 
be seen as a giant arc across the Galactic plane if it is pro- 
jected onto the Galactic coordinate system. The arc-like ap- 
pearance is due to the disky distribution of gas and stars 
inclined to the Galactic plane. The simulated arc-like belt 
can cover the observed OB stellar associations that are con- 
sidered to be formed in the GB, which suggests that the 
present model can reproduce well the distribution of young 
stellar populations in the GB. The DMC is currently located 
at {X, Y, Z) = (—6.3, —1.0, 1.4) kpc and the distances of 
the DMC from the Galactic center and the Sun are 6.5 and 
2.8 kpc, respectively. Fig. 3 shows that the projected distri- 
bution of older stars with stellar ages (fage) ranging from 25 
Myr to 37 Myr (corresponding to stars formed between r=8 
Myr and 20 Myr, i.e., during a starburst) with respect to the 
mass center of all stars appears to be ringed and moderately 



elongated. The distribution of the very young stars with 1 
Myr ^ fagc ^ 10 Myr also appears to be ringed, though 
a significant fraction of the stars exist inside the ring-like 
structure. The elongated morphology and size of the GB 
in the present simulation are similar to those in the pervi- 
ous dynamical models for the observed properties of the GB 
(Perrot & Grenier 2003). 

Fig. 4 shows that almost all stars have positive radial 
velocities (ur) with a mean of 7.1 km s~^, which means 
that the simulated stellar structure is expanding with re- 
spect to its mass center. The expansion rate of the structure 
at Ay — 400 pc is 17.9 km s"'^ kpc~^, which is quite similar 
to that (20 km s~^ kpc~^) derived from one of kinematical 
models of the GB (Lindblad 2000). Although there is a weak 
tendency that the outer stars are more likely to have larger 
v-c, the distributions of the stars in this l^Y-v^ phase space 
is highly inhomogeneous owing to random motion of stars. 
Fig. 4 also shows that stars with positive /S.Y are more likely 
to have positive Ux, which means that the simulated struc- 
ture has internal rotation with the direction consistent with 
that of the Galactic disk. The velocity gradient of for 
Ay| ^ 400 pc is 21.1 km s~^ kpc~^, which is comparable 
to the angular velocity of 7—24 km s~^ kpc~^ derived from a 
kinematical model based on observational data sets of stars 
in the GB (Lindblad 2000). 
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4 DISCUSSION AND CONCLUSIONS 

By using the standard formula for the timescale of coUi- 
sion/merging, the time scale of a cloud-DMC collision event 
(itm) can be estimated as follows (MaJiino & Hut 1997); 



where rih, u, and v are the mean rmmbcr density of the 
DMCs within the Galaxy halo, the geometrical cross section 
of gas clouds and a relative velocity between a DMC and 
a cloud. We here estimate rih for the central 50 kpc of the 
Galaxy (corresponding to the pericenter of the LMC orbit) 
and assume that v is the same as ono-dimcnsional velocity 
dispersion (= Vc/ \/(2), where Vc is the circular velocity thus 
220 km s"'^) of the Galaxy halo. If we use the results of the 
latest cosmological simulations, we can estimate the total 
number of DMCs (A*'h) with masses larger than ~ W^Mq 
within 50 kpc (Diemand et al. 2008). The total number of 
subhalos in the latest high-resolution cosmological simula- 
tion (Springel et al. 2008a) is by a factor of ~ 8 larger than 
that of Via Lactea II (Diemand et al. 2008). Considering this 
resolution effect, we adopt A^h = 100 rather than A'^ « 20 
(Diemand et al. 2008) in the estimation of tm- 

We here estimate tm for a DMC to collide with one of 
many GMCs (not with a specific GMC). Therefore a should 
be A^ci X TT X where A'ci and _Rci are the total number 
of the massive GMCs in the Galaxy and their sizes, respec- 
tively. Given that the Galaxy has ~ 4000 massive GMCs 
(Solomon et al. 1987) with masses more than 10^ M©, we 
can derive tm for iVci = 4000 as follows; 

*- = °-^^(M)"(life)"(l56]^)"«^^ 

This suggests that the cloud-DMC collision is not rare and 
thus that the dark impact can inevitably influence evolution 
of the Galactic disk to some extent. The collisions between 
DMCs with lower masses (Mdm ~ lO^M©) and less massive 
GMCs (with masses less than lO^M©) can be much more 
frequent than those described above, because a much larger 
number of these less massive DMCs are predicted to exist 
in the Galactic halo (Diemand et al. 2008). Those collisions 
between less massive DMCs and GMCs can form signifi- 
cantly smaller (thus less clearly visible) stellar structures in 
comparison with the GB. 

The present study has first demonstrated that the 
dark impact can dramatically change spatial distributions 
of forming stars and star formation histories within giant 
molecular clouds. This new mode of star formation is in 
a striking contrast to those resulting from small-scale tur- 
bulence within clouds (e.g., Larson 1981) and from global 
galaxy-scale dynamical instability (Elmcgrccn 1995) Future 
observational studies of stars in the GB by GAIA (Ferryman 
et al. 2001) will determine the three-dimensional structural 
and kinomatical properties of the GB stars and the individ- 
ual stellar ages and thus enable us to compare the results 
with those of the present study. This comparison will deter- 
mine whether and when the majority of stars in the GB were 
formed from a single coUisional event of a DMC in a more 
quantitative way. The GB is not just a magnificent stripe of 
stars in the night sky but also a fossil record containing pro- 
found physical meanings of a dramatic event which would 
have happened in the solar neighborhood. 
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